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Abstract 

We investigate whether dual strings could be solutions of the magnetohydrodynamics equa- 
tions in the limit of infinite conductivity. We find that the induction equation is satisfied, and we 
discuss the Navier-Stokes equation (without viscosity) with the Lorentz force included. We ar- 
gue that the dual string equations (with a non-universal maximum velocity) should describe the 
large scale motion of narrow magnetic flux tubes, because of a large reparametrization (gauge) 
invariance of the magnetic and electric string fields. It is shown that the energy-momentum 
tensor for the dual string can be reinterpreted as an energy-momentum tensor for magnetohy- 
drodynamics, provided certain conditions are satisfied. We also give a brief discussion of the 
case when magnetic monopoles are included, and indicate how this can lead to a non-relativistic 
" electrohydrodynamics" picture of confinement. 
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A conducting gas/fluid interacting with a magnetic field is described by magnetohydrodynamics 
(see ref. [lj for a review of this subject). A turbulent plasma has a strong tendency to become 
extremely intermittent, because the vorticity concentrates in thin vortex tubes, and the magnetic 
field concentrates in thin flux tubes. When one looks at computer simulations [g] of such a plasma, 
the situation looks very "stringy", in the sense of an intermittent network of "string"-like objects. 
This suggested to the author that it would be an interesting question to ask if dual strings (see e.g. 
the reviews mentioned in ref. |J) could be solutions to the magnetohydrodynamics equations. 

In the following we shall assume that the conductivity is infinite, and the viscosity vanishes. 
In computer simulations one sees [2] that the flux tubes are formed as parallel magnetic field lines 
over a long scale, but then the field lines spread out and the magnetic field becomes diffuse and 
weak. After some distance the field lines may again join in a flux tube. The diffuse behavior is 
not "string" -like, and is therefore from the very beginning not covered by our considerations. We 
do therefore not claim that all features considered in the following are realistic, but we hope that 
some of our results are of interest. 

In order to have a "stringy" situation, we consider flux tubes which are closed (or have an 
infinite length). We also assume that they move with the velocity of the fluid at any point. This 
puts severe restrictions on the dynamics to which our picture applies. In particular, we expect that 
the magnetic energy is roughly of the same order as the kinetic energy of the gas/fluid. 

One case where our considerations may be of direct use, is the case where somehow closed 
magnetic flux tubes are generated in the early universe, e.g. as defects in some phase transition. 
Then these flux tubes are expected to follow the Nambu-Goto string equations ||-[f|], at least to the 
first approximation. On the other hand, these tubes should also follow the magnetohydrodynamics 
equations, since such flux tubes interact with the charged particles. Thus, if these closed flux tubes 
are to keep their identity for some time (ultimately they are expected to decay), the Nambu-Goto 
string equations and the magnetohydrodynamics equations should be satisfied. Now the results in 
the following indicate that it is indeed possible to satisfy both set of equations at the same time. 

We start by considering a magnetic field in the form of an infinitely narrow flux tube 

B i (x,t) = bJ da^^ 5 3 (x-z(a,t)) . (1) 

The flux tube follows a curve given by Zi(a,t), where a is some parameter. The constant b has 
dimension ^/energy x length. Since we do not assume that the flux tube has a topological origin, 
b is not assumed to have a universal significance. We assume that the string moves with the same 
velocity as the fluid velocity v(x, t), i.e. 

«i(x, t) 5 3 (x - z(a, t)) = 6 3 (x - z(a, t)), (2) 

and B is assumed to be perpendicular to v. It is also assumed that the curve described by z(a,t) 
has no self-crossings or touchings, so that the tangent is unique. 

The non-universal quantity b has the interpretation of the magnetic flux through a surface 
perpendicular to B, as is seen by integrating the expression (|l]) over such a surface, 

J d 2 x ± B(x,t) = bj d 2 x ± J dz 5 3 (x - z) = b, (3) 
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where d 2 x± is perpendicular to B(x, t), and hence is also perpendicular to dz. The constancy of b 
is an expression of the well-known magnetohydro dynamics conservation of a magnetic flux which 
follows the motion of the fluid. 

Evidently eq. (|l|) represents the large scale behaviour of the magnetic field. It should be 
mentioned that an identical expression is used in hydrodynamics and in superfluid helium for the 
large scale behaviour of vortex lines (see ref. |J for a recent review). Here the vorticity u is given 
by 

w(x, t) = curl v(x, t) = const. / da ^ ^ - 5 3 (x — z{a, t)), (4) 

so the resulting velocity is of the Biot-Savart type. 
Taking the time derivative of eq. (JlJ) we get 

dBi(x,t) f d 2 Zj s3/ f dzidz k d 3 

b / d(J T^Z S (x-z) + b daj-^- j— 5 (x - z) 



dt J dtda J da dt dz^ 

u f j d2z i X 3, \ d(v k Bi) 

Using that 

(Vx(vxB)), = ^-^, (6) 
ox k dx k 



as well as eq. we obtain the equation 

dB 



0j V x (v x B), (7) 



where we used 



d(viB k ) d f dzidz k 3 f dzi d 3 f d 2 Zj 3 

" h l^r / ^-5T^r S (x-z) = -b da— —5 {x - z) = b / da—- <5 (x - z). (8) 



3xfc dxk J dt da J dt da J dadt 

Eq. (|7|) is the well known hydromagnetic form of the induction equation in the limit of infinite 
conductivity. Thus the "string" (i.e. narrow flux tube) (|l|) satisfies this equation for any curve 
Zi(a,t). 

The field B should also be divergenceless. It is quite easy to check that this is the case: 

f dz- d r d 

divB(x,t) = b / da—^ —S 3 (x-z) = -b / da— 5 3 (x-z) =0. (9) 
J da dxi J da 

Here we have used that the "string" is assumed to be closed, so that the contributions from the 
boundaries exactly cancel in eq. (§). 

In the limit of infinite conductivity the only existing electric field is the one induced by B, 
namely c E = — v x B, where c is the velocity of light, which cancels out in eq. (]?]). Now one of the 
main points of this paper is that the electric and magnetic fields can be unified in the tensor field 

W = - - ( 10 , 

Here r is an arbitrary "time" parameter, x^ in the delta function is the four vector (xq = t,x), and 
Zn is the four vector (zq,z). In the special parametrization ZQ^a, r) =r, we have 

Foi = —Bi, Fij = c eijkEk, (11) 
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Eqs. (0) and @ can also be unified in the single equation 



dx, 



0. 



fl2) 



Here we used the boundary condition that everything vanishes at r = +00 and r = — 00. The 
crucial point is now that the expression (|To|) is reparametrization invariant, i.e. F„ v is the same no 
matter which set of parameters (a, r) we use: A reparametrization, where a and r are replaced by 
a(a,r) and f(a, r), respectively, leaves invariant, and hence the physics does not depend on 
the parameters we use. The parametrization zq{&, t) = t used in the beginning is thus very special. 
In the general case z^(a, r) describes a surface in space and time, to be identified with the "world 
sheet" of the string associated with the (infinitely) narrow flux tube. If the time t is replaced by an 
arbitrary "time" parameter r, we thus have a very large "gauge" freedom, consisting of the freedom 
of choosing parameters a and r describing the same surface. Under very general conditions we can 
restrict this freedom by choosing an orthonormal system (see ref. ||). Now we have already taken 
the velocity perpendicular to the flux tube, i.e. 



dzi dzi 
~dt ~da 



0. 



(13) 



In order to have orthonormality we should further normalize the velocity and the tangent vectors 
according to 



Oz 



V 



+ 
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dz\ 2 
dt) 



1. 
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Here vq is a velocity introduced for dimensional reasons. In string theory the corresponding velocity 
is universal and equal to the velocity of light. Here we do not assume that vq is in any way universal. 
From eq. (|l~4| ) it follows that vq is the maximum (transverse) velocity the flux tube considered can 
have. 

It is now easy to see that the motion of the string is fixed to be given by the usual harmonic 
equation of motion for £j(cr, t). To see this, differentiate the normalization condition ( |i~4| ) with 
respect to a, 

d 2 Zi dzi 2 d 2 Zidzi 

+ v 01^jF— = 0- (15) 



dtda dt 



da 2 da 



From the orthogonality condition (13) we obtain by differentiation with respect to t 



d Zi dzi dzi d 2 Zi 



dt 2 da dt dtda 



0. 



Using this in eq. (^) we get 



which has the solution 



,d 2 Zj dz,: d 2 z-: dzj 



da 2 da dt 2 da 



(16) 



(17) 



d 2 2 



da 2 



Zi . / 2 • v dzi 



(18) 

where a is an arbitrary function. In the following we shall take a = 0. The function Zi{a,t) thus 
satisfies the harmonic equation together with the orthonormality conditions (13) and (|H|). The 
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solution to these non-linear equations is known from strings. The harmonic equation for z% has the 
solution ^ 

Zi(a,t) = qi+pit + i ^2 ~ exp \ inv °j) \ an ' i cos n J + ^ n >» sin n j) ' ( 19 ) 

n=— oo \ / \ / 

where I is some lenght scale, which does not have to be universal. The non-linear orthonormality 
conditions are solved by the conditions 

L n {a) = L n (f3) = 0, (20) 

where 

L n (x) = ^2 x n - m x m , x = a or (3. (21) 

Here the sum over m runs from — oo to +oo,and a summation over the space indices is understood. 
Also, we have defined Qo,i = /?o,i as beeing proportional to the constant pj. We refer to the literature 
for further discussion of the solution Q. 

We thus have the situation that the non-linear eqs. of motion can be solved in principle, and 
the motion of the string can be followed. In practice, however, the solution given in eqs. (19)-(^l|) 



is somewhat formal, and it should presumably be replaced by numerical solutions, where one starts 
from an initial curve and study the subsequent development. 

We have seen that the fundamental magnetohydrodynamics equation (]7|) is satisfied by the 
string ansatz (11). Now one of the central points of this paper is that if we make the ansatz 
([!]), then the existence of a more general underlying physical theory (Maxwell's) implies that the 
dynamics is really described by the field tensor (|l0[). This implies a very large freedom of 



parametrization, which can be restricted to be orthonormal (see ref. [3]) like in eqs. (13) and (14). 
Thus, any velocity field v(x, t), which is consistent with the ansatz (|l|) or (|i~0|), must necessarily 
satisfy the harmonic eq. (|l~8|). Thus it follows that either a pressure and a density field can be found 



such that the corresponding solution of the the Navier-Stokes equation also has the properties (13), 
@), and or the ansatz |§) or ^) is not consistent with magnetohydrodynamics. Encouraged 
by the fact that the ansatz (||) or (|i~0|) satisfies the induction eq. (0) , and by the result of computer 
calculations Q, we shall tentatively assume that the first possibility is realized. We therefore 
consider the Navier-Stokes equation, or rather the Euler equation including the Lorentz-force (since 
there is no viscosity term /iV 2 v), i.e. 

+ s,^ _ _sp + y x _ a i + B a 

Ot OXk OXi OXi 2 OXk 

Here -P(x, t) is the local pressure, and p(x, t) is the local mass density. For simplicity, we consider 
the non-relativistic version of the Navier-Stokes equation. Since Zi(a,t) in principle is known as a 
solution of the string eqs. (|K]), (|i~4"|), and ( |T8| ) with a = 0, this amounts to finding out the conditions 
imposed by equation ( |22|) on the pressure P and the density p. This is a non-trivial problem, since 
the gas/liquid may or may not penetrate the flux tube. We shall discuss this problem later. 

It should be noticed that the Navier-Stokes equation couples the magnetic field and the velocity 
dynamically, so the field B reacts back on the velocity field. Thus our previous statement that the 
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field B is supposed to move with the fluid in the induction equation (0) is to be understood in the 
sense that v is a functional of B. 

At this stage the reader may wonder why one cannot use the arguments leading to the Nambu- 
Goto equations @, @, and @ in the case of the vorticity (||). The difference is that e.g. in 
the hydrodynamics case the fundamental field is the velocity, and a; is a quantity derived from 
v. There is no physical point in considering a "unified" tensor (u,v x a;), and claim that this is 
somehow the fundamental tensor in hydrodynamics. This is in contrast to the Maxwell equations, 
which unifies E and B physically. Therefore, it has a physical meaning to exploit the large gauge 
invariance in F^ u , which is due to the reparametrization invariance, in much the same way that 
ordinary gauge invariance is often used to choose a simple gauge where the calculations are done 
in the most transparent way. 

We shall now compare with the Navier-Stokes eq. (|22|), which can be derived from the momen- 
tum flux density tensor (see Landau and Lifshitz Q, p. 215) 

IT fc = pviv k - B t B k + (P + -B 2 )S ik , (23) 

with 

at dx k ' y ' 

These equations are of course non-relativistic. We would like to compare Iljfc with the energy 
momentum tensor for the string. In order to do this, we need first to discuss the quadratic terms 
in B in eq. (23). In computing the quadratic terms from the ansatz (|l]) we encounter divergencies, 
due to the product of two delta functions. We therefore need a regulator, which we can take to be 
Gaussian, 

S 3 (x) « -4-7 exp , (25) 

where l± is the transverse dimension of a flux tube replacing the string, and where l± is assumed to 
be small relative to all other distances, of course. For simplicity we have regulated the longitudinal 
direction with the same regulator as the transverse directions. Regulators of the type p5[ ) have 
been used in the study of superfluid helium ||. 

Let us now consider one of the quadratic terms on the right hand side of eq. ( f23|) . Formally we 
have 

B k Bt = b 2 J dada'z k {<T,t)z' i (a',t) 5 3 (x - z(a,t))5 3 (x - z(a',t)), (26) 

where the prime on z denotes the derivative with respect to the first argument. If we use this 
expression, we have an infinity. If we insert the regulator (|25|), we can use that the integrals have 
support only when x = z(cr, t) = z(<r', t), to within the distance l±. Thus, to the first approximation, 
we can replace one of the delta functions by its value at zero argument. If o"o is the parameter 
which satisfies x = z(o"o,i), then the requirement that x ~ z(a,t) to an accuracy of order l±, is 
satisfied for \gq — a\ ~ l± <C I, as is seen by expanding the solution ( |i~9| ) around a®. Thus the 
effective range of the a- integration is of order l±. Up to a numerical constant we thus obtain 

b 2 f 

B k Bi « -j— / daz' k (a, t)4(a, t) 5 3 (x - z{a, t)). (27) 

7T 2 / | J 
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The moral lesson of this exercise is that we cannot accept the ansatz (Q) as an exact string equation, 
but we must always think of the string as a narrow flux tube. 

Now we shall show that the magnetohydrodynamics equations based on the ansatz (|l|) and 
discussed above, can be obtained by a reinterpretation of the string energy-momentum tensor. The 
latter has been used often in the study of cosmic strings (see ref. Q), and is given by 

^ ring = -t- I dadr - z'A 5\x - z(a, r)), (28) 

where we use four-vector notation, and where the dot indicate the derivative with respect to r. 
The constant in front is introduced for reasons to be made clear in the following. Like in eq. (|10|), 
T^ rmg is a priori defined in terms of the arbitrary parameters a and r, but in the expression ( p8| ) 
we have already used the orthonormal " gauge" for simplicity (see the string literature Q and Q for 
a discussion of this point). It is easy to see from the harmonic eq. (18) that there is conservation 
of energy and momentum flow, 

^string 

^ = 0. (29) 
dx^ 

The boundary terms disappear because the string is closed and because we assume that everything 



disapppear at infinite times. Now from eq. (28) we have, by doing the r-integral 

Tt n& = I da^(x - z(a, t)) = pv u (30) 

where we have defined 2 

p&) = -^hr~, I da6 ^ x ~ z( - a > '))■ ( 31 ) 

It is easy to see that p satisfies the continuity equation, and that it has the dimension of a mass 
density. We also have 

b 2 f 

T kf mg = P v k y i TV / daz 'k z 'i # 3 (z - z(<?, *)) = P v kVi - B k Bi. (32) 

7T2 lj_ J 

Here we used eq. (27) in the last step. It is now seen that the magnetohydrodynamics equation in 
the form (24) is equivalent to the string energy- momentum conservation ( |29| ) with v = i, provided 
the density of the fluid inside the string is identified as in eq. ( |3l|) , and provided we take 



P + ^B 2 = const. (33) 



This condition means that the space-space components of two tensors ( |23| ) and (pq) differ by a 
constant only. Such a difference is without dynamical importance. 

It is interesting that the condition (^) is the same as the variation of the pressure in the exact 
hydromagnetodynamic wave discussed by Landau and Lifshitz (Q-see p. 221). Eq. (|33|) tells us 
that the hydrodynamic pressure should have a dip inside the string (flux tube). 

Assuming the ansatz defined by eqs. (|T|) , p|) , (|3l|) , and (|33|), we therefore have derived the 
result that the dual string energy-mometum tensor can be reinterpreted as the energy- momentum 
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tensor of magnetohydrodynamics in the limit of infinite conductivity. The classical motion of the 
dual string can thus be understood as a motion in a fluid under the influence of a magnetic field 
generated by the string itself as in eq. (JlQj) . In this connection it should be remembered that the 
string cannot be considered as an ideal zero-width string, but must be given a small transverse 
dimension. Of course, the same is true for cosmic strings, where eq. ( p8| ) is also used 01, where the 
global properties are well described by this string energy-momentum tensor. 
A simple, but not very interesting, closed string solution is given by 

z(<7, f) = I cos ujt (cos— ,sin— ), (34) 



I ' I 

where I is the lenght scale introduced previously, and where the maximum velocity is given by 
vq = lu. For < cot < 7r /2, this solution represents a circular flux tube with decreasing radius, 
corresponding to a fluid which moves towards the origin. It is more interesting that many results 
are known about dual strings (see refs. 0-0], where further references can be found), so we hope 
that the re-interpretation mentioned above can be of use in magnetohydrodynamics. Conversely, 
in cosmology it may be of interest that magnetic strings follow the magnetohydrodynamics way of 
life. 

We end this discussion by emphasizing that string-like (or flux tube) solutions have previously 
been considered by many authors, as is clear from the literature quoted in ref. O. A recent review 
by Semenov is given in ||. He also derives harmonic equations, but they differ from ours, since in 
his case the term z" is multiplied by the density p, and orthonormality is not used. His solution is 
more general than ours, since the condition (|3^) is not imposed. In general, the energy-momentum 
tensor corresponding to Semenov's equations cannot be reinterpreted in terms of dual strings. This 
shows that there are other string-like solutions of the magnetohydrodynamics equations than the 
dual type considered in the present paper. 

So far we have considered closed strings. It is, however, quite easy to add monopoles. The 
induction equation ([?]) is then changed to 

dB 

~dt 

where the magnetic current is given by 



V X (v X B) - jmagnetic, (35) 



jmagnetic = b z(l,t) 5 3 (x - z(l, t)) - b z(0, t) 5 3 {x - z(0, t)) • (36) 

Here we inserted the ansatz (|l|). The divergence equation (||) is of course also changed, 

div B(x, t) = -b 5 3 {x - z(l, t)) +b5(x- z(0, t)). (37) 

The constant b can be fixed by Dirac's quantization condition. Thus the generalised magnetohydro- 
dynamics equations have magnetic strings between the magnetic monopoles, which may therefore 
be confined. Of course, we do not know if there are non-confining solutions which are more favoured 
energetically. For the sake of the argument, let us assume that these strings are favoured. Then 
we have a pedagogically simple picture of quark confinement: We know that Maxwell's equations 
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are invariant under dual transformations. Hence, we can imagine a situation where we have a 
"universe" (vacuum) filled with magnetic monopoles, with a "dual Ohm's law", 

j magnetic — ^magnetic B, (38) 

where tT magne tic is the conductivity of the monopoles, and where eq. (|38|) is valid in the rest frame. 
Assuming now that the magnetic conductivity is infinite, we would have "electrohydrodynamics" 
which is dual to the usual magnetohydrodynamics. Thus the B-field vanishes in the rest frame of 



the fluid. Thus, for example, instead of eq. (35) we would have 

— = V X (v X E) -jelectric, (39) 

at 

since the only magnetic field is given by c B = v x E in the moving frame. Eq. d3^) is therefore 



one of the Maxwell equations. The electric current has an expression analogous to eq. (36). It 
should also be assumed that the "dual displacement current" dH/dt vanishes. The string solution 
of eq. ( |39| ) is thus an electric string, 

E = b' /dcr^M 5 3 (x-z(a,t)), (40) 



da 

where b' is a constant. Here the field z(cr, t) follow the ortonormal dual string equations. There 
is an equation analogous to eq. ( |37| ) for div E. The electric charge, which is related to b' , is thus 
confined. 

This scenario of quark confinement is somewhat similar to the dual superconductor picture Q. 
It is perhaps easier to imagine electrohydrodynamics than it is to imagine a dual superconductor, 
because the former is a simple consequence of the duality inherent in the Maxwell equations. On 
the other hand, it is not clear if our simple picture can be generalized relativistically, and hence it 
can only be regarded as some low energy approximation. In any case, a crucial ingredient in both 
cases is that conductivity should be very large and monopoles are needed in the vacuum. In QCD it 
may be possible to distinguish the two cases in lattice gauge calculations, e.g. by checking whether 
the dual London equations are satisfied ||. In electrohydrodynamics these dual superconductor 
equations are not expected to hold. 

I thank Ake Nordlund for an illustrative demonstration of his amazing computer simulations of 
magnetic flux tubes, and I thank Axel Brandenburg for communicating his knowledge and insight 
in magnetohydrodynamics in an almost infinite number of conversations. Thanks are also due to 
Mark Hindmarsh for pointing out the existence of ref. ||. 
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